The bound on viscosity and the generahzed second law of thermodynamics 
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We describe a new paradox for ideal fluids. It arises in the accretion of an ideal fluid onto a black 
hole, where, under suitable boundary conditions, the flow can violate the generalized second law 
of thermodynamics. The paradox indicates that there is in fact a lower bound to the correlation 
length of any real fluid, the value of which is determined by the thermodynamic properties of that 
fluid. We observe that the universal bound on entropy, itself suggested by the generalized second 
law, puts a lower bound on the correlation length of any fluid in terms of its speciflc entropy. With 
the help of a new, efficient estimate for the viscosity of liquids, we argue that this also means that 
viscosity is bounded from below in a way reminiscent of the conjectured Kovtun-Son-Starinets lower 
bound on the ratio of viscosity to entropy density. We conclude that much light may be shed on 
the Kovtun-Son-Starinets bound by suitable arguments based on the generalized second law. 

PACS numbers: 51.20.-(-d,66.20.-f d,04.70.Dy,11.15.-q,11.25.Hf 



I. INTRODUCTION 

The use of general principles to investigate systems 
whose microscopic makeup is unclear can be very reward- 
ing. Sometimes this method gives information on a whole 
class of systems. Among such principles, thermodynam- 
ics, believed to be of universal applicability, stands out. 
An instructive example of its use is the application of the 
second law of thermodynamics to the problem of viscous 
flow, for which it permits the inference that the two vis- 
cosity coefhcients must be positive without need to resort 
to microscopic expressions for the latter [l| . 

In this work we use the generalization of the second 
law of thermodynamics holding in the presence of black 
holes (the generalized second law — GSL) to reach fur- 
ther conclusions about the shear viscosity coefficient of 
an arbitrary fluid. We do this by describing a new 
paradox for ideal fluid flow in the presence of a black 
hole. This indicates that the correlation length of a 
real fluid cannot be arbitrarily small. By implication 
the energy-momentum tensors describing systems which 
display macroscopic fluid behavior must be subject to a 
restriction: the shear viscosity, a function of the thermo- 
dynamic state of the system, cannot be arbitrarily small. 
Thereby the GSL opens an alternative macroscopic ap- 
proach to the recently proposed lower bound on viscos- 
ity 

The GSL is a unique law of physics that bridges ther- 
modynamics and gravity. It is rooted in the understand- 
ing that a black hole, basically a puregravity entity, is 
endowed with well deflned entropy ^, proportional 
to its surface area. The GSL |5|, |8j] then claims that 
the sum of the entropy of all black holes and the to- 
tal ordinary entropy in the black holes' exterior never 
decreases. While this formulation is reminiscent of the 
ordinary second law, the GSL is exceptional in that it re- 
lates ordinary entropy — a rather elusive object from the 
mechanical viewpoint — and the surface area of the black 
hole (formally the area of its horizon) whose evolution is 
quite mechanical in nature. From this point of view it 



is little surprise that the GSL has provided unexpected 
information on entropy. 

An example is the upper bound on the entropy of 
weakly self-gravitating thermodynamic systems (the uni- 
versal bound on entropy — UBE [9|,[l3|)- While in partic- 
ular cases the bound can be verifled directly, it is the 
GSL which really makes understanding of the bound in 
generic situations easy. However, the GSL gives more 
than just a simpler way to see some results derivable by 
other means. At the microscopic level the GSL repre- 
sents a piece of the yet to be established theory of quan- 
tum gravity. In particular, this law permits, in principle, 
to draw conclusions that from the microscopic viewpoint 
would only be derivable from a fundamental theory com- 
bining quantum mechanics with gravity. For example, 
the GSL gives an indication of the number of particle 
species in nature [13, [llj. In another example, the GSL 
reveals an a priori bound on the strength of the electro- 
magnetic interaction (l^ . 

Let us now describe the ideal fluid paradox revealed 
by invoking the GSL. We assume the existence of physi- 
cal fluids with arbitrarily small shear viscosity 77 at fixed 
values of some two thermodynamic variables, say the en- 
tropy and the energy densities, s and p, respectively. (In 
this work we consider only simple fluids, so the values of s 
and p completely determine the thermodynamic state of 
the system.) Here "physical" means, among other things, 
that the fluid satisfies the GSL; this could, in principle, 
be checked by considering the fluid at the microscopic 
level, but we shall not go into such detail. 

The assumption that the fluid has arbitrarily small vis- 
cosity allows us to describe its flow as ideal fluid flow, 
possibly with shocks (the zero viscosity limit of a given 
flow can be non-trivial; see Refs. [l|,[ll). It turns out, as 
we shall show, that for sufficiently slow accretion of the 
fluid onto the black hole the overall entropy decreases 
and the GSL is violated. The realizability, in principle, 
of such slow accretion flows will be demonstrated, so that 
the small viscosity assumption engenders a paradox. 

Of course, ideal fluid paradoxa exist already in non- 
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relativistic physics, for example, the famous d'Alembert 
paradox. This maintains that an ideal fluid with no 
boundaries exerts no force whatsoever on a body moving 
through it with constant velocity. In particular, there is 
no lift force, so swimming or flying would be impossible 
in such a fluid. Thus the established fact that swimming 
is possible in any fluid implies that every real fluid must 
have nonzero viscosity. There is however an essential dif- 
ference between the d'Alembert paradox and the paradox 
described in this work. D'Alembert's paradox says noth- 
ing about the actual magnitude of 77. In contrast, our 
paradox implies, as we shall see, the existence of a lower 
bound on 77 for given s and p. 

This conclusion is in concord with the recent ingenious 
conjecture of Kovtun, Starinets and Son (KSS) [Ij^ (see 
also Ref. [4j]). Based on holographic calculations of the 
viscosity coefficient for certain strongly coupled quantum 
field theories with gravity duals, they suggested that the 
viscosity rj of a general, possibly non-relativistic, fluid is 
subject to the miiversal restriction 
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Currently this bound is considered a conjecture well sup- 
ported for a certain class of field theories — see the de- 
tailed discussion in Ref. and the references therein. 
In the sequel we discuss the relation between the ideal 
fluid paradox presented in the next section and the KSS 
bound. We also argue that a frequent objection to the 
validity of the KSS bound is likely to be ruled out by the 
GSL. 

Our paper is structured as follows. In Sec.[TT]we show 
explicitly that the slow accretion of a truly ideal fluid 
onto a Schwarzschild black hole leads to a contradiction 
with the GSL. One escapes from the paradox by recog- 
nizing that every fluid must have a nonvanishing corre- 
lation length which restricts the range of applicability of 
the ideal fluid paradigm. In Sec. IIIII we obtain a lower 
bound on the correlation length and a generic estimate 
of the viscosity of real fluids, which together bound the 
viscosity to entropy density ratio from below. Although 
this is not yet the KSS bound, we consider there the con- 
nection between it and the UBE, and argue that the GSL 
provides a natural frame for elucidation of the origin of 
the former. In the Sec. lIVI we summarize our results and 
arguments. The realizability, in principle, of the slow ac- 
cretion flow assumed in Sec. |lT] is demonstrated in the 
Appendix. 

Unless otherwise stated, we work in units with c = 
h = k = 1 , where c is speed of light and k is Boltzmann's 
constant. Our metric signature is (—,+,+,+). 



II. IDEAL FLUID PARADOX 

In the present section we consider ideal fluid flow into 
a spherical black hole. For some flows, we demonstrate 
that it is possible for the GSL to be violated so that the 



total entropy of the system decreases. It follows that the 
assumption of a perfect continuum down to an arbitrarily 
small scale is not consistent with the GSL. 



A. Entropy balance in accretion onto a black hole 

Consider a flow, not necessarily spherically symmetric, 
in which fluid is absorbed by the black hole. The rate of 
change of the total entropy S of the system is the sum 
of the rate of change of the entropy of the black hole 
exterior, S^xt, and that of the black hole entropy Sh- 



dS_ 



dt 



dt 



(2) 



Here and below we use Schwarzschild coordinates 
{x^ , x"^ , , x^) = (t,r,9,4>), with r > rn where rn is 
the Schwarzschild radius. 

We first calculate dSext/dt. Let the fluid's proper en- 
tropy density be s. Since the fluid is assumed ideal, there 
is no dissipative contribution to the entropy current den- 
sity which is thus purely convective, and must take the 
form sU^, where t/^ is the fluid four- velocity. The fluid 
can carry entropy into the hole leading to a decrease of 
Sext- The explicit expression for this comes from the 
entropy balance equation P, [H, [l6| , 



gsU^) = 0, 



(3) 



where g 



' sin'^ 9 stands for the determinant of the 



Schwarzschild metric g^i^. 

dSext/dt of the black hole exterior {r > th) is thus 



d_ 

dt 



I sU°y^drd9d(P = - drd9d(pdr{sU''y/^) 



(4) 



We have not included a contribution from the outer 
boundary of the domain of integration because we intend 
to specialize to stationary flows. In any such situation the 
entropy flow into the hole per unit t-time is given by the 
r.h.s. of Eq. ([4]). The expression is non-positive because 
< for infalling matter. We have assumed that the 
flow is differentiable, and that it contains no shocks; oth- 
erwise there exists an additional contribution to dSext / dt 
associated with the entropy generation in the shocks [l| . 

Let us now consider the second term of the r.h.s. in 
Eq. The absorption of the fluid by the black hole 
increases the latter's mass Af , producing an increase in 
the black hole entropy given through energy conservation 
by 



dS 



H 



H- 



dt 



dM 
~dt 



dE, 



dt 



(5) 



where is the black hole temperature, and M is its 
mass. We shall now write down the flux of energy into the 
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black hole. Let = {1, 0, 0, 0) be the Killing vector asso- 
ciated with the stationarity of the black hole. Then the 
energy-momentum tensor of the fluid, T^'', must obey 

iCT^l;^ = 0, (6) 

or equivalently 

d.iV^T,^) = 0. (7) 
Since T^^—g is minus the energy density of the fluid, 



'H- 



dt 

drded(f)dr{T^- 



dt 



rjir 

-t n 



-gdedcf). (8) 



by 



The energy-momentum tensor of the ideal fluid is given 



(9) 



where p is the energy density and p is the pressure in 
the comoving frame. Consider now the normalization 



condition J/'^C/^ = — 1, written as 



(10) 



Since ^Jgee and ^Jg^V^ are physical velocity com- 
ponents, they should be bounded at r = r}j. Hence 
since is future and inwardly pointed, and grr ^ oo as 
r vh, 'we can infer from the last equation that Uq = 
at the horizon . Combining all the above we find 



at Ih Jr=rH 



(11) 



Thus dSn /dt > in harmony with Hawking's area theo- 
rem 

We observe from Eqs. (jl]) and (fTT|) that while the (neg- 
ative) rate of change of Sext is proportional to the first 
power of U"^, the (positive) rate of change of the black 
hole entropy is proportional to the second power of U^. 
Thus for sufficiently small U^{rH), the total entropy of 
the system will decrease, in violation of the GSL. Explic- 
itly we have 



dS 
dt 



\U'\{P + P)' 



Th 



(12) 



where we stress that at the horizon is never positive. 
We observe that when the accretion velocity Uac (the 
suitable mean value of \U^\ over the horizon) obeys 



p + p' 



(13) 



the total entropy decreases and the GSL is broken. We 
now proceed to search for such flows. We first reconsider 
known explicit solutions. 



B. Validity of the GSL for the Bondi flow 

The above formulae apply for a generic, not necessar- 
ily spherically symmetric, accretion flow onto the black 
hole; in particular it may have non- vanishing and W^. 
As an example, consider the well known Bondi flow (see 
Ref. [3 and references therein). Bondi flow starts from 
rest at infinity and is spherically symmetric. For a poly- 
tropic equation of state, 



p — Kn , 



(14) 



where V is the adiabatic exponent, the accretion velocity 
is given by Ugc « 1 for F 7^ 5/3 and Uac ~ 0.782 for 
F = 5/3 is"]. It follows that Bondi flow wiU obey the 
GSL if the following condition holds: 



< 



p + p 



UirrH, F ^ 5/3, 

[(0.782) X 47rrff, F = 5/3. 



(15) 



We have used the usual expression Th = (47rr//)^^ for 
the Schwarzschild black hole. 

The above conditions would seem to obtain for any 
fluid with positive pressure thanks to the UBE 9]. This 
bound states that the entropy Sb of a generic, weakly 
self-gravitating, thermodynamic system satisfies 



1^ < 2nRb, 



(16) 



where Eb is the body's total energy, while i?f, is its linear 
size. Our experience is that the inequality here is usually 
a strong one. 

Now the minimal size of a parcel of fluid consistent 
with the fluid description is its correlation length I. For 
a gas, I is the mean free path, while for a liquid it is 
typically the intermolecular distance. Applied to such a 
parcel the UBE tells us that (see the next Section for 
more details) 



s/ p < 2ttI, 



(17) 



where we have passed from the total entropy and energy 
to their densities by dividing by the parcel's proper vol- 
ume. Again, this bound will in most cases be a strong 
inequality. Since the continuum description down to ac- 
cretion at the black hole makes sense only it I <^ rn, we 
may conclude that the inequalities p3|) will always hold 
with at least one order of magnitude difference between 
the r.h.s. and the l.h.s. of the equations. In the above 
argument we have tacitly assumed that p > 0. Negative 
pressure is often discussed in cosmology (dark energy); 
the functioning of the GSL in the face of dark energy is 
fraught with subtleties p^ . 

Thus Bondi flow satisfies the GSL by virtue of the 
UBE. Clearly, this happens because the accretion veloc- 
ity approaches the speed of light, Uac ~ 1- We now turn 
to examples of slow ideal fluid accretion flows that do 
violate the GSL. 
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Paradox in slow accretion flow of ideal fluid 



relativistic version of the Bernoulli equation reads [ij 



It is clear from the previous analysis that the GSL 
holds for any accretion flow of ideal fluid with Uac ^ 1- 
We expect such flows whenever spherical accretion starts 
from a distance from the black hole large compared 
to (for Bondi flow = cxd). But what is Uac in the 
limit ri th'^ We now show that Uac vanishes in that 
limit, thereby giving an example of ideal fluid flow that 
violates the GSL. 

We study steady, spherically symmetric accretion flow 
which starts at r — ri with the initial flow and sound 
velocities, Ui and a^, specified there. We assume that 
the relevant initial conditions are physically realizable; 
see the Appendix for the discussion of this assumption. 
Following Ref. [ll| we introduce u = -U"^ > 0. The 
continuity and the Euler equations can be written as 



D 



D ' 



(18) 



where a prime denotes 9/9r, n is the baryon number 
density, and 



D 



u2 - (1 - 2M/r + u'^)a'^ 



Di = - 

n 



(l^2M/r + u^)—-- 



D2 



2u^/r - M/r^ 



(19) 
(20) 

(21) 



with a the local sound velocity. Obviously D{r = th) = 
{uln){l — a^) is always positive due to the causality con- 
straint < 1 (sound velocity is smaller than velocity of 
light) ^Ts*!. Here we consider only flows with Uac > so 
u 7^ at the horizon. 

In Bondi flow D becomes negative at large r which sig- 
nifies that there exists the so-called sonic radius such 
that D{r — r^) = 0. As is clear from Eqs. (US]), where D 
stands in the denominator, the sonic radius is a special 
though regular point of the Bondi flow, which in other 
situations could signify the presence a shock. By con- 
trast, here we shall choose a range of initial parameters 
for the flow which ensure that D{r = ri) > 0, so that 
there is no sonic point. We assume initial parameters for 
the flow satisfying 



2M 



< 1, a- < 1. 



2M\ , 
1 U? 



, 2M 
^ < < 1 , 

n J Ti 



(22) 

(23) 
(24) 



It is easy to sec from Eq. that D{ri) is indeed posi- 
tive. 

We do not need to derive the explicit form of the flow 
in order to flnd Uac- It suffices to consider the conserva- 
tion of the Bernoulli integral along the streamlines. The 



P + P 



1 - 



2M 



u ~ const. 



(25) 



Assuming for simplicity the polytropic equation of 
state one finds the following relation 



(26) 



r - 1 - a2 



where m is the baryon mass. The above relation allows 
us to rewrite the Bernoulli equation as 



1 - 



2M 



= const. 



(27) 



Evaluating the above equation first at r — and then 
at r = r;, equating the results, using inequalities (j22p - 
(j23p and assuming that at the horizon remains much 
smaller than unity (which we have verified numerically), 
we find that 



Ua 



1 - 



2M 

n 



(28) 



where the ellipsis stands for subleading terms. The above 
expression holds for any K and F, the latter assumed to 
be not too close to unity. 

Actually, the use of the polytropic equation is not es- 
sential. Using Eq. ((25| directly we would find 



Ua 



njrH) [p{ri)+p{n)] 
n{ri) [pirn) +p{rH)] 



1 



2M 



(29) 



If we may assume that there are no abrupt changes in 
{p+p)/n throughout the flow, we recover a result of form 
Eq. (pS)) . We may conclude that an ideal fluid accretion 
flow which starts close to the horizon with initial pa- 
rameters complying with inequalities (|22p - (|23p can have 
arbitrarily small Uac- According to the criterion it 
will thus incur a violation of the GSL. 

Of course for the above argument to be convincing, 
we must still demonstrate that it is possible for a flow 
to start near the horizon with sufficiently small veloc- 
ity (and sound speed). In the Appendix we show in 
detail that a cord which respects fundamental physical 
constraints can be used to bring objects to rest arbitrar- 
ily near the horizon. One can envision the setting up of 
the initial conditions for the fluid we require at r = ri by 
this means. 



D. How to remove the paradox? 

The above violation of the GSL raises a paradox. In 
the real world the GSL cannot be violated. This law has 
been shown to follow from fundamental concepts in quan- 
tum theory and classical gravity [2(|. Can the paradox 
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be removed without calling on new physics just for the 
occasion? 

An obvious solution would be to call on Hawking's ra- 
diance to generate entropy that would at least compen- 
sate for the decrement of total entropy pointed out in 
Sec. Ill CI After all, in free Hawking emission the ther- 
mal radiation entropy creation rate exceeds the associ- 
ated rate of decrease of the black hole entropy [2l| . But 
there exist reasons to reject this as the resolution of our 
problem. The Hawking radiation is capable of preventing 
a violation of the GSL that would arise if high entropy 
radiation of the same nature, i.e., electromagnetic, were 
injected into a black hole with Th above the effective ra- 
diation temperature [2l[. Its efficacy here is related to 
the principle of detailed balance in equilibrium (in fact, 
for incoming thermal radiation at temperature Th the 
Hawking radiation exactly balances the entropy decre- 
ment). Now detailed balance refers to modes of the same 
physical system. There is no detailed balance between 
radiation and fluid modes. Thus if entropic radiation 
were injected alongside the fluid accretion, the Hawking 
radiation might prove incapable of compensating for the 
entropy decreases of the two kinds. 

The above discussion also suggests we should look for 
a way out of the paradox that hinges on the physics of 
the fluid itself. In our discussion in Sec. Ill Cl we assume 
that one can deposit the fluid at an arbitrarily small dis- 
tance from the horizon at an arbitrary velocity. Since 
the physical paradigm used is that of fluids, one should, 
as a matter of principle, demand that the said distance 
is still larger than the correlation length I at which the 
hydrodynamic description first becomes applicable. 

To be in the Schwarzschild spacetime a small proper 
length I <C M away from the horizon at th = 2M corre- 
sponds to the Schwarzschild coordinate r — 2M + MSi ; 
here i5; <C 1 and is given explicitly by 



Si 



8M2 



(30) 



Let us thus substitute 2M -f M 5i for rt in our expres- 
sion 



1 - 



2M 



I 

4M 



(31) 



Thus the correlation length limits the smallness of the 
accretion velocity Uac- 

With Eq. ([3T|) the expression in square brackets in 
Eq. takes the form 



27Tl{p+p)] 



(32) 



where we used the expression for the black hole tempera- 
ture. In order for the GSL to be obeyed the above factor 
must be nonpositive. This is actually guaranteed by the 
fluid version of the UBE, Eq. pT]) , as long as the pressure 
p is nonnegative (see the next Section for more details). 
This removes the paradox. The moral of the discussion 
is that a paradox arises if one relies on the continuum 



description of a fluid down to an arbitrarily small scale, 
that is if one takes the notion of ideal fluid literally. To 
be rid of paradoxa one must take cognizance of the flnite 
correlation length of any physical fluid, and must accept 
that the entropy capacity of fluid matter is limited ac- 
cording to bound p7|). 



III. GSL AS FRAMEWORK FOR 
ELUCIDATING THE KSS BOUND 



The paradox uncovered in Sec. Ill CI and its resolution 
in Sec. IIIDI clearlv show that the ideal fluid paradigm is 
inconsistent. In particular, the picture of a fluid as a per- 
fect continuum is shown to be physically unacceptable: 
the fluid in question must have a nonvanishing correlation 
length. The medium is a fluid only over scales exceeding 
the correlation length. 

It turns out that flniteness of the correlation length 
is incompatible with the vanishing of various transport 
coefficients like shear viscosity and heat conductivity, an- 
other feature of the ideal fluid paradigm. For a gas there 
is a simple way to see this. The usual estimates of the 
mentioned transport coefficients [23j have them propor- 
tional to the mean free path of the gas' molecules. But in 
a gas the mean free path and correlation length are the 
same thing. Thus the finite correlation length forces the 
transport coefficients of the gas, in particular the shear 
viscosity, to be nonvanishing. We shall see in the sequel 
that the same conclusion applies generally to any liquid 
as well. Thus a fluid can be fully compatible with the 
GSL only if it is dissipative to some extent (and thus not 
ideal). Another way of putting this is that a fluid with 
arbitrarily small shear viscosity is unphysical. 



A. Lower bound on correlation length 

How large must the correlation length I be? Micro- 
scopically speaking it must obviously be at least as large 
as the intermolecular distance. But can we say something 
without delving into the structure of the fluid? Let us ap- 
proach this question in the spirit of Wilson's work on the 
renormalization group (2^ . We consider some thermody- 
namic system with the typical linear size R. Because the 
system is macroscopic (by definition), one can represent 
any extensive thermodynamic variable as some quantity 
times the system's volume V. For example, one can write 
the entropy as S" = sV and the energy a.s E — pV. Now 
decrease R. For sufficiently small R, already comparable 
with the system correlation length Z, the system ceases 
to be macroscopic, and extensivity is generally lost. This 
means that the expressions for the entropy and energy 
densities themselves become dependent on R. The sys- 
tem is no longer a continuum. Can one set a lower bound 
on the correlation radius I purely from macroscopic con- 
siderations? 

The standard answer to the above question would be 
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in the negative. But, surprisingly, the true answer is yes. 
As we have seen, for a macroscopic system the UBE can 
be restated in the form (fl7|) , which immediately leads to 
the inequality 



log r] 



I > 



2Trp 



(33) 



The main point here is that a macroscopic system with 
size much smaller than sjl'Kp would violate the UBE; 
thus the correlation length cannot be much smaller. We 
emphasize, again, that I may easily be much larger than 
the minimal scale ([33| . The above result answers the fol- 
lowing question: given a macroscopic system with given 
entropy and energy densities, what is its minimal possible 
correlation length? It shows that macroscopic quantities 
do "know" about the minimal correlation length of the 
system. 



B. Lower bound on shear viscosity 

We mentioned earlier that a nonvanishing correlation 
length I ensures that the transport coefficients of a gas 
(which are proportional to /) do not vanish. For exam- 
ple, for shear viscosity one has the order of magnitude 
estimate [231 



77 ~ p/a, 



(34) 



where the speed of sound a is of the order of molecular 
speed |23]. The above estimate together with Eq. 
give that Ty/s is subject to the inequality 



(35) 



Let us show that the estimate and thus the 

bound l|35)) . must hold also for liquids sufficiently far 
from a critical point. In this preliminary treatment we 
neglect heat conductivity and bulk viscosity. We note 
that, for any fluid, density perturbations at scale much 
larger than I generate sound waves. On the other hand, 
perturbations with scale much smaller than I are not co- 
herent and produce no sound. Furthermore, one can still 
use hydrodynamics asymptotically to describe the evolu- 
tion of perturbations with scale I. Then the demand that 
there is no well-defined sound at smaller scales gives the 
asymptotic condition that at scale I the wave decay time 
~ pl^/r], as found from hydrodynamics fH, should be 
comparable with the wave period ^ l/a. This produces 
the estimate (IMl) for liquids. 

Alternatively, the results ([M)l - (P5|) can be recovered by 
considering a sound wave already propagating through 
the liquid. Its Fourier components with reduced wave- 
lengths near or below I should decay over a distance com- 
parable to I since we cannot have macroscopic flow at 
those smaller scales. Now from the macroscopic point of 
view, the decay must be caused by transport processes 
controlled by the viscosities or heat conductivity. Ac- 
cording to fluid theory [J, a sound wave of wavelength 
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FIG. 1: Plot of logio of experimental viscosities of eleven 
pure liquids (data taken mostly from Ref. [23) vs logio of the 
estimates from Eq. (I34p with I identified with the average in- 
termolecular separation. Both viscosities and estimates are in 
SI units (mPa s) . Nonstandard symbols are "M" for methane, 
"P" for propane, "A" for acetone, "E" for ethanol, "U" for 
undecane (C11H24), "Ni" for nitrobenzene (pheniI-N02) and 
"G" for glycerol (CaHgOa). A repeated symbol correspond to 
viscosities at several pressures and temperatures. The solid 
line is the locus of = pal; the dotted lines demarcate the 
region where the viscosity lies within a factor of 3 of pal. 



A penetrates a distance of order pa[\/2iTY /rj into the 
liquid before damping out. According to the above argu- 
ment, for (A/27r) < I this penetration length should be /, 
which gives Eq. ([M)) again. 

For relativistic fluids the above consideration should 
be modifled slightly. Here the wave decay time includes 
p -\- p instead of p , which leads to 



rj ^ {p + p)la 



(36) 



For most realistic fluids p < p/3, see [25j, and p + p ^ 
p, so this modiflcation is not essential for the order of 
magnitude estimate of 77, see the next subsection for an 
example. 

The above estimates, however, must fail sufficiently 
close to a critical point. Both j] and I diverge at the crit- 
ical temperature with the power law 77 cx holding in 
the vicinity of the critical point. Here A is much less than 
1 [13] . This is incompatible with Eq. ([M]) because both p 
and a are finite at the critical point. The reason for the 
failure of the estimate is complications in the asymptotic 
matching procedure above, related to the difference of 
the critical exponents of the various quantities involved. 
Another issue is that near a critical point the static cor- 
relation length I becomes different from the scale beyond 
which hydrodynamics applies. 

Nevertheless, sufficiently far from the critical point the 
estimate ([M)) seems to produce remarkably good results, 
as shown by Fig. [T] for a number of pure liquids. (The 
only resounding failure is glycerol, one of the most vis- 
cous fiuids known. The problem with it may be that / is 
much larger than the intermolecular separation, which we 
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routinely used as estimator for I.) This fact is remarkable 
because it allows to reach conclusions on the magnitude 
of viscosity — quite possibly of practical value — even for 
liquids for which no microscopic theory is available. 

Let us now turn to the lower bound on the viscosity- 
entropy density ratio ([35]) . It evidently holds if the esti- 
mate works. Note that Eq. ([55)1 will most often be 
a strong inequality, particularly for nonrelativistic sys- 
tems, for which the UBE is known to be a very liberal 
bound. Thus the bound (|35p is generally valid far from a 
critical point. On the other hand, near the critical point 
the bound ([55)1 holds trivially because its l.h.s. diverges 
while its r.h.s. remains finite. Thus the divergence of 
shear viscosity at the critical point [111 only strengthens 
the bound. 

The considerations in this subsection provide a strong 
argument in favor of the existence of a generic lower 
bound on the viscosity-entropy density ratio. The re- 
formulation of the arguments with regard to the other 
transport coefficients, including thermal conductivity, is 
left for future work. 



C. The quark-gluon plasma 

It was the challenge presented by the quark-gluon 
plasma (QGP) which motivated the activity leading to 
the formulation of the KSS bound. As a practical use 
of our estimate let us apply it to the QGP. Quan- 
tum chromodynamics (QCD) predicts the existence of 
such a form of matter where the usually confined quarks 
and gluons are essentially free forming a fireball with a 
typical size ~ 10 fm. The transition from the hadronic 
state of matter to the QGP occurs at densities of 1 — 10 
GeV/fm'^. QCD lattice calculations predict that the 
transition should take place at a critical temperature 
Tc ~ 160 MeV, and that the speed of sound for a QGP 
above T « 300 MeV is about a - 0.3 c (sec Rcf. [28] and 
further references therein). 

Measurements on the QGP at Brookhaven's relativistic 
heavy ion collider (RHIC) indicate that the ratio of en- 
ergy density e to entropy density s is roughly proportional 
to the QGP temperature T (see, e.g.. Fig. 3 in Ref. [s^p. 
Specifically, for the high energy results of RHIC one has 



0.8 kT. 



Using this and estimate (|34p we have 



k 7] 
s h 



keal 0.8 kTal 



0.2 



(37) 



(38) 



where in the last step we assumed a ^ 0.3 c at kT — 500 
MeV, / ~ 0.3 fm corresponding to a number density n ~ 
40 /fm^, and employed he ~ 200 MeV fm. 

Our estimate for the QGP viscosity is thus in harmony 
with the KSS bound ([1]). Moreover, it is in quite rea- 
sonable agreement with the experimental data which re- 
quire a shear viscosity to entropy density ratio as low as 

rj/s < 0.2 Hi, [13, mi, m. 



D. KSS bound — a tightened entropy bound for 
fluids 



We draw attention to the similarity between inequal- 
ity psp with c and h restored, namely 



2tt 



(39) 



and the conjectured KSS bound ([T]). At first sight this 
inequality seems to fall well below the KSS bound since 
in many cases a <^ c. However, we must keep in mind 
that, especially in such nonrelativistic circumstances, we 
expect the inequality to be a strong one. The minimum 
rj/s may thus be well above the literal bound p9|) . and 
not necessarily at variance with Eq. ([T]). For relativistic 
media a ~ c, and there is no difference to speak of be- 
tween bounds ([!]) and ([39| . Thus arguments based on 
the UBE, and ultimately on the GSL, seem to suggest a 
raison d'etre for the mysterious KSS bound. 

Before proceeding we wish to provide an alternative 
viewpoint. The KSS bound Eq. ([T]) is usually interpreted 
as saying that a fluid cannot be too perfect. However, in 
the form 



s < 47r?7, 



(40) 



the bound is really an entropy bound, specifically an up- 
per bound on the entropy density of an arbitrary fluid (we 
return to the use of units with c = ?i = 1 as in the previ- 
ous sections). Clearly the above bound is reminiscent of 
the UBE, Eq. and its fluids version, Eq. P7)) . 

But at least for non-relativistic fluids, the KSS bound 
is a tighter entropy bound by orders of magnitude than 
the UBE. Restoring the speed of light c we introduce the 
variable 770 = Icp (l being, again, the correlation length, 
either the mean free path for a gas, or, typically, the 
intermolecular separation for a liquid) which has the di- 
mensions of viscosity. Thus Eq. pT|) gives s < 271770- Us- 
ing the estimate ([M)) we have 77/770 ~ a/c. It follows that 
for nonrelativistic fluids Eq. (|40p is indeed a much tighter 
entropy bound than Eq. ([T7]). Thus the KSS bound can 
be viewed as a tightened version of the UBE in guise of 
Eq. (fT7)) which is available for the class of systems ex- 
hibiting macroscopic fluid behavior. 

Both the ideal fluid paradox of Sec. Ill CI and the in- 
terpretation of the KSS bound as an entropy bound for 
fluids suggest the GSL as a natural frame for investigat- 
ing relations between entropy bounds. Because the UBE 
is obtained in the most transparent way via the GSL, one 
may hope that the KSS bound should likewise be obtain- 
able most simply with help of the GSL. We turn now to 
considerations in this direction. 



E. Where does the KSS bound come from? 

First we remark that it is not clear how to obtain the 
KSS bound directly from microscopic physics. Inspec- 
tion of the Green-Kubo formula 34] for the viscosity in 
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terms of fluctuations shows no apparent connection of the 
viscosity and the entropy. Such consideration affords no 
special status to the ratio ry/s, at least not from the view- 
point of nongravitational physics. On the other hand, 
the KSS conjecture emerged from holographic type argu- 
ments that connect quantum field theory with gravity. It 
might thus turn out that a derivation of the KSS bound 
requires use of the still nonexistent theory of quantum 
gravity. But even if this were true, one need not loose 
heart. There is general agreement that black hole en- 
tropy, and the GSL which hinges on it, reflect some as- 
pect of quantum gravity. 

The KSS is not universally accepted. Many objections 
to it rely on scenarios where very many particle species 
are supposed to exist in nature (see Ref. [14| for an exam- 
ple and references). These are the same objections raised 
against the validity of the UBE (see Ref. for refer- 
ences). In fact, any entropy bound whatsoever invites 
attacks of this sort, for if the number of species that may 
show up is unlimited, the entropy can be made as large 
as desired while keeping parameters like energy or total 
particle number fixed. Although it is formally true that 
with many species available, the KSS bound must fail, 
this does not detract from its heuristic usefulness; one is 
often interested on the entropy of a specific system, one 
whose particle content is fixed ahead of time. It should 
also be mentioned that many particle scenarios, if not 
arbitrarily legislating particle proliferation, conjure up 
fine-tuned or baroque setups to beget the required large 
number of species of quasiparticles or excitations. This 
aspect considerably decreases the appeal of the species 
proliferation arguments. 

Both of the above considerations make it clear that a 
ground up approach to deriving the KSS bound is un- 
likely to succeed. Alternative, indirect approaches are 
needed. A point in favor of employing the GSL to in- 
vestigate the origin of the KSS bound is that the GSL 
"knows" the actual number of species in nature [13, [lH . 
For example, black hole entropy, which plays a crucial 
role in the GSL, should in principle depend on the num- 
ber of elementary fields, yet all derivations endow it with 
a fixed coefficient which may be thought as determined 
by the actual number of species. This feature would act 
to neutralize the above mentioned argument against the 
validity of the KSS bounds. 



IV. SUMMARY AND OUTLOOK 

We have worked out the expression for the rate of 
change of the total entropy of a system consisting of a 
Schwarzschild black hole and ideal fluid which can ac- 
crete onto the hole. The well known Bondi flow is the 
case of fast accretion with velocity near to the speed of 
light; for it the GSL is always obeyed due to the UBE. 
For sufficiently slow accretion velocity the flow violates 
the GSL. The question of whether flow with the required 
small accretion velocity is a realistic option is answered 



in the affirmative by a critical revision, in the Appendix, 
of the venerable argument [35| claiming that it is impos- 
sible to adiabatically lower mass to near the black hole 
horizon. 

Since it is known from microscopic considerations (20| 
that any physical system should comply with the GSL, 
the above result would constitute a serious paradox if the 
ideal fluid is a continuum, as usually considered. We find, 
however, that the paradox can be defused if one takes into 
account that the continuum picture must break down at 
some level by virtue of the fluid having a nonvanishing 
correlation length. An auxiliary role in the nullification of 
the paradox is played by the universal entropy bound. It 
gives a lower bound on the correlation length in relation 
to the fluid's entropy density. 

The lower bound on the correlation length of an ar- 
bitrary thermodynamic system in terms of its (macro- 
scopic) entropy and energy densities is an unexpected and 
thought-provoking aspect of the UBE. It makes it clear 
that the GSL, though a macroscopic law, "knows" about 
the microscopic structure of matter. The lower bound 
must rank as one of the most impressive consequences 
of the extension of the second law of thermodynamics to 
include black holes. 

The breakdown of the continuum description of a fluid 
has momentous consequences: discreteness of matter and 
thermal fluctuations necessarily engender nonideal be- 
havior parameterized by the viscosity and heat conduc- 
tion coefficients. We have shown that one can expect 
the existence of a universal lower bound on the viscosity 
of an arbitrary fluid. Our argument yields a seemingly 
novel estimate for the viscosity of a fiuid far from the 
critical region; this estimate is shown to be better than 
an order of magnitude estimate for a range of pure liq- 
uids. In addition, the estimate works reasonably well for 
the quark-gluon plasma, which viscosity was a subject of 
much study lately [H [sl IM HI . 

Together with the lower bound on the correlation 
length, the viscosity estimate yields a lower bound on vis- 
cosity per unit entropy density which may be expressed in 
terms of the sound velocity of the fluid. The appearance 
of the speed of sound here, in contrast to the (implicit) 
speed of light in the original KSS arguments, may signify 
that the KSS argument can be improved, or may perhaps 
disclose that the field theories considered by KSS as rep- 
resentatives of strongly coupled fiuids are not sufficiently 
generic — in all these theories the sound velocity is com- 
parable with the speed of light. An additional issue is 
whether the bound can be extended to other transport 
coefficients. 

If there is indeed a lower bound on the ratio of viscosity 
to entropy density for an arbitrary fluid, then it is likely 
to be rooted in some very basic principle of physics. We 
have shown that the generalized second law of thermo- 
dynamics, combining as it does gravity with thermody- 
namics, can be such a principle. Starting from it, one can 
anticipate that the viscosity to entropy ratio of an arbi- 
trary fluid obeys a certain lower bound. In particular, the 
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GSL incorporates physics that weakens the recent objec- 
tions tM] to the bound's very existence. On the basis of 
the above, we propose that the use of the GSL may lead 
to a clarification of the KSS bound's origin. Although we 
have made no concrete progress towards reaching a final 
sharp inequality, we have given for the first time physi- 
cal arguments indicating that this elusive bound, whose 
existence is sometimes disputed, indeed exists. 
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APPENDIX A: EQUILIBRIUM OF MATTER 
NEAR A BLACK HOLE AND REALIZABILITY 
OF SLOW ACCRETION FLOWS 

Gibbons [11] claimed that one cannot lower a rope to 
near a black hole at an arbitrarily small velocity, without 
it being torn at some finite distance from the horizon. 
The conclusion was based upon the assertion that the 
stresses that arise in the rope become so strong that they 
violate the weak energy condition [36l . [37| which must 
be obeyed by the energy-momentum tensor of physical 
matter. Here we reexamine the problem and show that, 
in fact, the energy condition need not be violated and 
thus that adiabatic lowering down to the horizon is, in 
principle, possible. 

In the stationary state of some matter in the vicinity 
of a spherical black hole, the energy momentum tensor 
obeys 



1 5(t;, 



1 dg^, 



2 ^x^^ 



(Al) 



where we used the symmetry of T^^ 15]. Below we con- 
centrate on configurations for which Tj^ is diagonal. Us- 
ing the same coordinates as in Sec. Ill A| the /i = com- 
ponent of the above equations is satisfied automatically, 
while the radial component gives 



1 9(r;r^) 



MTI 



r^{l-2M/r) r^{l-2M/r) 



0. 



(A2) 



The above equation gives a linear relation between the 
different components of the stress tensor which need not, 
a priori, respect the weak energy condition. As an ex- 
ample, consider a thin spherical shell in equilibrium. For 
the shell the contribution of in the above equation 
is negligible. This can be seen by noting that van- 
ishes at the shell boundaries, and as a result its values 



within the shell vanish together with the ratio of shell 
thickness to radius. Far from the hole, where the shell 
is describable by the classical linear elasticity theory, the 
maximal value of is found to be proportional to the 
square of the above ratio. So the thin shell is basically 
supported by the tangential stresses Tg and T^, which 
must be equal because of spherical symmetry. Eq. (|A2p 
gives \T^\/\T^\ = M/(2r - 4M). Now the weak energy 
condition would demand ]Tg ]/]Tq] < 1. It is thus clear 
that a stationary thin shell would violate the condition 
at r < 5Af/2. Thus a physical thin shell, i.e. one obey- 
ing the energy condition, cannot support itself arbitrarily 
close to the horizon (see Ref. [sl] for the more detailed 
discussion). 

Ref. [35| argued that the rope cannot be in equilibrium 
with its lower end arbitrarily close to the black hole, sim- 
ilarly to the thin shell above. Here we wish to correct 
this conclusion. We shall assume that the rope fibers 
can be considered radial, such as in the case of a coni- 
cal rope filling the portion of space defined by some solid 
angle. Since the rope can only support stresses along its 
fibers, the general form of its energy-momentum tensor 
is T't = diag[-p, S, 0, 0], and Eq. gives 



d [r'^S) _ M(p + S) 
dr ^ 1 - 2i\f/r 



(A3) 



The above linear equation can be used to express the 
stress S{r) in terms of the rope density p and the bound- 
ary condition at the lower end of the rope, r = Tq, which 
is defined by the load. The solution can be written as 
a sum S{r) = Si{r) + S2{r), where Si{r) describes the 
stresses caused by the load, while S2{r) describes the 
stresses caused by the rope's own weight. The expression 
for Si (r) (the solution for a weightless rope) expresses the 
"constancy of the tension" along the rope: 



Siir) 



rg^l-2Af/ro 
- 2Af/r 



S{ro), 



(A4) 



where >5'(ro) is determined by the weight of the load; for a 
point mass m we have r§-\/l — 2M/ro S{ro) = — m. Note 
that Si{r) is negative, as befitting a tension. The "con- 
stancy of the tension" gives monotonically decreasing 
]S'i(r)] as one moves out because the rope cross-section 
increases with r by the assumed symmetry of the fibers. 

The expression for 5*2 (r) (the solution without the 
load) describes how the force caused by the rope above 
some r balances the weight of the rope below: 



S2{r) = - 



1 



^ , dr'^m:^. (A5) 

For definiteness we shall consider below the case of the 
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rope with a constant density po where integration gives 



52 (r) - - 



r Vl - 2M/r 



(A6) 



Vl-2M/ro + Mln ^ , 



We now show that the above sohition does not violate 
the weak energy condition |S'|/po < 1, and in particu- 
lar, that no violation occurs even when the rope's end is 
arbitrarily close to the horizon. By choosing sufficiently 
large po we may always disregard Si{r), so it is sufficient 
to show that max |5'2(r)|//9o < 1- It is easy to see from 
Eq. (jASP that the maximum of |<S'2(r')| = — S'2(r), at given 
Po and M, grows as ro decreases. This just means that 
the maximal stress is the bigger the closer the lower end 
of the rope is from the horizon. Therefore, it is enough 
to show that max | >5'2 (r) | /po < 1 for tq which is infinitesi- 
mally close to th — 2M. In general, 15*2(^)1, vanishing as 
it is at both endpoints, r = tq and r = oo, has a unique 
maximum at r» e (ro,oo). It is easy to see numerically 
that in the limit vq rn the point also tends to rn- 
Defining e = (1 — 2Af/r*)^/'^ we find analytically 



2M _ 8e^ 
ro 3 
max 152(01 p. 1 _ 

PO 



(A7) 



follows from the above that the weak energy condition is 
obeyed for any rg > r^/. Thus, at least from the view- 
point of the energy conditions, it is possible in principle 
to lower a body adiabatically all the way down to the 
horizon by means of a suitably constructed rope . 

Our result is at variance with that of Gibbons [ssl ]. 
The discrepancy arises because he missed a term when 
calculating a certain 4-divergence (Eq. (5) in Ref. [ssj). 
By taking this term into account, the main result, Eq. (8) 
in Ref. [35] , now reads 



dT _ dV_ 
T + a ~ ~V 



(A8) 



Here T = SA is the tension, A the rope's cross section, 
a — pA the energy per unit proper length of the rope, 
and V = \/—Ka K°- is the norm of the timelike Killing 
field A'" (the spacetime is taken to be stationary). In the 
original derivation, the term T in the denominator of the 
l.h.s. of Eq. ((X8|l is absent. The corrected Eq. is 
indeed in accord with our Eq. (jA3[) . Observing that in 
the Schwarzschild geometry V = \/l — 2M/r and further 
A~^dA/dr = 2/r for the case of a conical rope, one easily 
checks that Eq. (fA8|) reduces to our Eq. ((A^ . 



where e ^ 1 and the dots stand for subleading terms. It 



Our analysis demonstrates that no problem of principle 
militates against the creation of the boundary conditions 
required for the slow accretion flow considered in Sec. Ill CI 
here. A suitably designed rope could be used to deposit 
parcels of fluid at rest near the horizon. 
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